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Abstract 

We investigate the index of the Neuberger's Dirac operator in abehan gauge theories on 
finite lattices by numerically analyzing the spectrum of the hermitian Wilson-Dirac operator 
for a continuous family of gauge fields connecting different topological sectors. By clarifying 
the characteristic structure of the spectrum leading to the index theorem we show that the 
index coincides to the topological charge for a wide class of gauge field configurations. We 
also argue that the index can be found exactly for some special but nontrivial configurations 
in two dimensions by directly analyzing the spectrum. 
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It is generally believed that the nontrivial topological gauge field configurations are responsi- 
ble for nonperturbative aspects such as the large r)-rj' mass splitting in QCD and the fermion 
number violation in the standard model. Lattice gauge theories are the most promising 
approach to such nonperturbative phenomena. However, the configuration space of lattice 
gauge fields is topologically trivial since any configuration can be continuously deformed into 
the trivial one. This is not quite disappointing because physically interesting gauge fields are 
smooth and we can impose a kind of smoothness conditions that restrict plaquette variables 
within a small neighborhood of unity The space of smooth link variables acquires a non- 
trivial topological structure and it is possible to define topological invariants geometrically 
in terms of lattice gauge fields |[|, ^ |]. The index theorem plays a role of a touchstone 
in lattice QCD |, |, 0, |. 

In the case of Ginsparg- Wilson (GW) Dirac systems |1^, [Tl|] it is possible to define 
exact chiral symmetry [|1^ and the index theorem on the lattice relating the index of the 
GW Dirac operator to chiral anomaly is known by Hasenfratz, Laliena and Niedermayer 
(HLN) [|T3], |T^. Their index theorem is applicable not only for any gauge field configurations 



but also for any GW Dirac operators. Since the topological structure of the gauge fields is 
refiected on the very construction of the GW Dirac operator rather implicitly, the HLN index 
theorem tells very few about the relationship between the index and the topological invariant 
of gauge fields for strictly finite lattices. Nevertheless we expect that a lattice extension 
of the index theorem relating the index directly with the topological invariants of gauge 
fields can be established for sufficiently smooth configurations and physically acceptable 



GW Dirac operators \T^\. In fact it is possible to relate the axial anomaly (the index 
density) to the topological charge density for abelian theories on the infinite lattice quite 
generally by invoking general principles such as the gauge invariance and the locality |T^. On 
strictly finite lattices, however, we cannot appeal to the locality. This makes the differential 
geometrical method developed in ref. inapplicable. 



In this note we focus our attention to the index of Neuberger's Dirac operator [|T^ for 
compact U(l) theories on finite periodic lattices in two and four dimensions and aim at 
establishing the index theorem for smooth background gauge fields.[| To achieve the goal we 
investigate the spectrum of the related hermitian Wilson-Dirac operators numerically for a 
family of link variables connecting configurations with distinct topological charges and clarify 
the characteristic structure of the spectrum leading to the index theorem. Such an analysis 
has already been reported in ref. within the framework of the overlap formalism. We 
shall extend the analysis from the point of view of Neuberger's Dirac operator to convincing 
extent. We then show some analytic results concerning the spectrum and the index. They 
are very useful to understand the numerical results. 

Let us begin with the definition of Neuberger's Dirac operator on a d = 2N dimensional 
euclidean hypercubic regular lattice L'^ with periodic boundary conditions. For simplicity 
we choose the lattice spacing a = 1 and take the lattice size L to be a positive integer. We 

'''The chiral anomahes have been estabUshed in ref. [Il6| in the continuum Umit. See also ref. Il7|. 
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first introduce the hermitian Wilson-Dirac operator H by 

HH^) = 7d+i|(rf - m)^{x) - X: [^^U,{x)ij{x + A) + ^^U,{x - fi)ij{x - A)^ | , (: 

where m is the fermion mass and the Wilson parameter r = 1 is assumed. The 7-matrices 
are taken to be hermitian and satisfy {7/x,7z/} = 25^,^. We use 7^+1 = (— z)2 7i ■ ■ - 7^. The 
link variables Ufj,{x) are subject to the periodic boundary conditions ?7^(a; + LO) = f/^(x). 
We also employ the periodic boundary conditions for ?/'(a;).0 The Neuberger's Dirac operator 
D is then given by 

D = l + 7d+i^. (2) 

The m is chosen in the range < m < 2 to avoid species doubling. We shall often use 
m = 1 in the following analysis and denote H{1) simply by H. Whenever the m dependence 
is concerned, we write it explicitly as H{m). 

For the Neuberger's Dirac operator the index theorem of HLN can be stated as 

indexD = Tr7,+i - ]^D^ = "^Tr-^ . (3) 

This equation relates the index of D to the spectral asymmetry of H since the trace on 
the rhs is nothing but the number of positive eigenvalues of H minus that of negative ones 
[p^ . It enables us to analyze the index of the complicated operator D in terms of the much 
simpler operator H that is ultra local on the lattice. 

On finite lattices D is well-defined unless det H = 0. By excising gauge field configura- 
tions leading to det H = the space of the link variables becomes disconnected. We shall 
refer to the decomposition of the configuration space as the analytic. The index theorem 
associates indexD to each connected component of the space of link variables. 

The lattice index theorem (j^) is valid for any gauge field configuration with det if 7^ 0. 
However, we can restrict ourselves to smooth configurations by imposing 

sup||l-P^,(a;)||<r/, (4) 

where P^^ is the standard plaquette variable and 77 is a positive constant to be mentioned 
later. The norm \ \A\ \ of a matrix A is defined by 

\\A\\= sup \\A^\\, (5) 

\M\='^ 

where \ \(f\ \ is the standard norm for vectors. It is well-known that the space of such smooth 
configurations still contains any physically interesting gauge field configurations with non- 
trivial topological structure if the size of the lattice is taken to be sufficiently large [2D|] . 

^We may take anti-periodic boundary condition in the d-th lattice coordinate. This does not lead to any 
significant change in the numerical analysis. 
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For suitable choices of t] the space of hnk variables becomes disconnected and one can 
define topological invariants geometrically in terms of the link variables. We shall refer to 
the decomposition of the space of link variables as the geometric. We describe this in detail 
for compact U(l) theories with which we shall be exclusively concerned. 

In the case of compact U(l) theories (^) can be cast into the simpler form: 



sup < r/' . (r/' = 2sin 

r TT for = 2 and D < r;' < 7r/3 for > 4 



(6) 

We may choose < r/' < tt for = 2 and < r^' < 7r/3 for > 4 g |T3|, ^ ||. The space of 
link variables is decomposed into a finite number of connected components characterized by 
a set of magnetic fluxes <^y,y = 2TTm^i^ through /iz/-plane defined by ^0 



L-l 

4>fii' = ^f^'^i^ + sfi + tv) , (7) 

s,t=0 

where m^i, = — m^^ is an integer and is the abelian field strength given by 

1 . 



Ff,u{x) = -lnP^^(x) . {\F^u{x) \ < n) 



Note that the rhs of (0) is independent of x by virtue of the Bianchi identity d[xF^i^] = 0. 
Here 9^ is the forward difference operator defined by d^f{x) = f{x + fi) — f{x). The (p^i, 
is not only a smooth function of the link variables but also is a topological invariant. This 
implies that any two gauge field configurations with distinct sets of magnetic fluxes cannot be 
deformed into each other continuously and, hence, belong to different connected components. 
Conversely, any gauge field configurations satisfying (Q) can be continuously deformed into 
one another if they have the same set of magnetic fluxes ||20|. This implies that the space of 
gauge fields with a given set of magnetic fluxes that satisfy (|]) are connected. In particular 
any gauge field with a given set of magnetic fluxes = 2T:m^y can be continuously deformed 
to a configuration with constant field strengths F^y{x) = (p^i^/L"^. For later convenience we 
denote the space of link variables that satisfy and give a set of magnetic fluxes (p^y by 



We can define a topological charge by the lattice analog of Chern character |T^, ^ as 



X • ■ • X F^j^^j^{x + /ii + z>i H h fXN-i + i>N~i) , (9) 

where is the Levi-Civita symbol in d dimensions. Qat is a smooth function of the link 

variables within a connected component and takes an integer value given by [0, ^] 
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A bound on \Qn\ can be easily found from (H) by noting (y) as 

We thus find two different definitions for decomposition of the space of hnk variables, the 
analytic and the geometric. The analytic decomposition is based on the requirement det H ^ 
and indexD is associated to each connected component. The geometric decomposition is 
based on the condition and Qn can be assigned similarly to each connected component. 
In general they are not mutually consistent and one can find gauge field configurations that 
satisfy (^) and det H = simultaneously. However, if we choose rj within the range 



H{= H{1)) cannot have zero-modes pT| , and any connected component of the space 
of link variables satisfying is contained in some connected component obtained by the 
analytic decomposition. Therefore it is very natural to expect that the index of D given by 
and the topological charge @ coincide with each other. The precise form of the index 
theorem for abelian gauge theories can be stated as : 

For the link variables satisfying ^ and p^ j the index of Neuberger's Dirac operator D and 
the topological charge Qjy are related by 

indexD = {-ifQN ■ (13) 

To establish (|l^) on an arbitrary connected component W^^,^ it suffices to show the equality 
for a particular gauge field configuration with constant field strengths F^y{x) = (p^y/L"^. 

One way to find indexD for an arbitrary configuration is to count the net number of 
spectral flows of H{m) as m varies from to 1 [0, 52]. Since H{m) is known to have a 



symmetric spectrum for m < [|T8|, the net number of eigenvalues that changes the sign 
from plus to minus just coincides to the index of D at m = 1 as can be easily seen from (j^). 

In what follows we shall take another path since we intend to understand the behaviors 
of indexD under continuous changes of link variables that include topological jumps. We 
consider a one-parameter family of link variables {Uj^^}o<t<T connecting two gauge field 
configurations {f/^°^} and {Uj^^}, where both topological charge and index are known for 
{Uj^^} but only topological charge is known for {Uj^^}. If {Uj^^} and {U^^} belong to two 
distinct connected components, Uj^^ goes through configurations violating det H 7^ or (|^). 
We expect jumps in indexD and Qat at some values of t. In the actual analysis we shall 
choose U^^ = 1. The indexD for {Uj^^} can then be found by counting the net number of 
the eigenvalues of H that change the sign from plus to minus as t varies from to T. 

In two dimensions we consider the one-parameter family of link variables given by 



Ui\x) = exp 



-tt—X2dxi,L-l 



1/2^ {x) = exp 



2tt 

itj^Xi 



(0 < < L - 1) (14) 
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The field strength is then computed as 

Fi'^i^) =^J2- - n)6,,,L-iSx,,L-i for < t < , (15) 

where n is an integer with \n\ < — 1 and is chosen to satisfy \Fi2{L — 1, L — 1)| < vr. 
The topological charge Qi is nothing but mu and is given by Qi = n. The hnk variable 
( p!4D connects uniform field strength configurations Fi2\x) = 27m/L^ of topological charge 
Qi = n {n = ±1, ±2, ■ ■ ■) continuously to the trivial one Fi2\x) = for which H has a 



symmetric spectrum. As a function of t, Qi has a discontinuity at t = = 
(j = 0,±l,±2,---). 



top _ (2j + 1)L' 
' 2(L2 - 1) 



Before going into detail about the numerical results, it is helpful to note the following 



facts: (1) The eigenvalues A of H{m) are bounded by |A| < d + \d — m| (2) The 

eigenvalue spectrum of H{m) is periodic in t. (3) The eigenvalue spectrum of H{m) is 
symmetric with respect to the point t = 0, i.e., —A is an eigenvalue of H{m)\t^_t if A is an 
eigenvalue of H{m). From (2) and (3) we see that indexD is an periodic odd function of 
t and vanishes at t = 0, ± jl, where the spectrum is symmetric. 

We have analyzed the eigenvalue spectrum of H. over the range —LP'/2 < t < for 
2 < L < 15 by searching for the values of A such that Re[det(i/- A)]Re[det(if-A- AA)] < 
for small enough AA.| The eigenvalues that are degenerate or approximately degenerate with 
an even multiplicity cannot be detected in this method. This, however, causes no problem 
since for most values of t the eigenvalues are not degenerate. 

In Figure |l| the whole spectrum of H is shown for L = 6. The characteristic structure of 
the spectrum is not changed with the lattice size L. It can be summarized as follows: 

• Most of the eigenvalues are contained in the upper and the lower trapezoid regions 
symmetrically separated by the parallelogram region. 

• Every time t increases by unity from an integer an eigenvalue belonging to the lower 
trapezoid crosses the parallelogram upward and moves to the upper trapezoid. In par- 
ticular H has zero-modes for a sequence of values of t on the interval of the horizontal 
axis cut by the parallelogram. They are almost equally spaced. 

• For integer values of t a kind of clustering of the spectrum occurs on the trapezoids and 
there appear large gaps on the parallelogram. In particular, there are 2sL eigenvalues 
with degeneracy r at t = r if L is a product of two positive integers r and s. 

Let K,{L) be the minimum value of t (> 0) where an eigenvalue changes the sign from plus to 
minus, then n{L) is larger than This can be qualitatively understood by noting that 

^In the actual computation the approximated values of dct{H — A) have imaginary parts. However, the 
ratios of the imaginary part to the real part are very small and the real part is considered to well approximate 
the determinant. 
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Figure 1: 

the smoothed boundary of the parallelogram is convex and goes through the points (0, 1) 
and (1/^/2, —1). As we shall mention, it is possible to show analytically that for even L the 
eigenvalues of if at t = satisfy 

V^-1<|A|<1, V5<|A|<v^+l. (16) 

Hence the point (L^/2,— 1) can be regarded as a lower right corner point of the upper 
trapezoid. We see that there is a large gap —y/b < A < — 1 on which H has no eigenvalue. It 
is analogous to the gap — 1 < A < 1 for the free theory and can be observed at integer t. As 
for the lower bound on n{L) we also note the exact result indexD = — for L a multiple 
of 4. They support the aforementioned observations based on the numerical analysis. 

On the interval —k,{L) < t < k{L), 2n{L) eigenvalues flow upward and change the sing 
at fj"^ (j = -n{L) + 1, -n{L) + 2, ■ ■ ■ , n{L) - 1, n{L)) from minus to plus, where tf'^ < tjf*^ 
for j < k. The indexD jumps by —1 at these points and monotonically decreases from the 
maximum n{L) to the minimum —n{L). Since n{L) k,{L) for large L, the maximal index 
that can be realized by grows faster than as L increases. In Table |l], n{L) is given 
for 2 < L < 15. We find numerically L"^ /n{L) ^ 3.6 for large L. 

Numerically t™'^ is approximately equal to t^"^ for |j| < n(L) and t™*^ > t^°^ for j > 0. 
In the case of L = 4, for instance, fj"^ for j = 1, ■ ■ ■ , 4 are given by 0.5387 (0.5333 • ■ ■), 
1.6203 (1.6), 2.7161 (2.666 ■ ■ ■), 3.8447 (3.7333 ■ ■ ■)• The difference tf'^-tf'^ is monotonically 
increasing in j for a fixed L, while it is monotonically decreasing in L{> 4) for a fixed j{> 0). 
We have found that the maximal deviations t]i^(L) ~ ^n°(L) scattered at around 0.2 and fj''^ 
for |j| <^ f^{L) approaches to as L becomes large. The numerical values of t™*^ — 
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Table 1: 



for fixed j can be well-fitted to a curve ajL ^ for large L, where the coefficient aj depends 
on j. In the limit L — > oo, both f^^ and approach to j ' + -. This is in accord with the 
observation that on the infinite lattice the topological invariants can be investigated quite 
generally by using the geometric concepts [|T3|. In Figure |^ index/? and —Qi are plotted 
over the region < t < L^/2 for L = 5. The correspondence between the index and the 
topological charge is very excellent on the region \t\ < k,{L). We thus conclude that indexD 
coincides to —Qi for gauge fields sufficiently close to configurations with the constant field 
strength |Fi2| < 27rn(L)/L^. This is consistent with the numerical results given in ref. . 

The jumps of indexD and those of —Qi, however, do not exactly synchronize and the 
equality between indexD and —Qi is violated near the discontinuities, where the field 
strength (^) at the corner point Xi = X2 = L — 1 takes a large value ~ ivr. The fi- 
nite deviations between t^"*^ and imply that the analytic decomposition of the space of 
link variables into connected components by imposing det H ^ and the geometric one by 
imposing (|^) are not consistent with each other for a general choice of r^. If r/ is chosen to 
satisfy 0<r7<l — 1/ as given by (|T2|) , the condition excludes uniformly the regions 
where the jumps both in the index and in the topological charge occur. We thus find that 
index/} coincides with —Qi for any gauge fields contained in and the index theorem (|13|) 
holds true in two dimensions. On the configuration space Uj^^tj, the maximal index is given 
by [r/'L^/27r]. In Table |I| we give bounds nf^{L) on the index for 2 < L < 15 for comparison 
with n{L). 

The spectral flow crossing A = is complicated on the interval n{L) < t < L^/2. The 
eigenvalues flow both downward and upward. The downward flows are dense and concen- 
trated in narrow regions, where the index increases rapidly. Then the index decreases slowly 
due to the sparse upward flows until the next dense downward flows approach to the hori- 
zontal axis. The downward flows and the upward flows are repeated alternately depending 
on the size of the lattice. Finally we have indexD = at t = I? jl. Similar thing also 
happens on the interval —L'^jl < t < —k,{L). We should consider, however, such peculiar 
behaviors on the outer regions k{L) < \t\ < L^/2 as a kind of lattice artifacts and the cor- 
responding field configurations as unphysical. In other words the field configurations should 
be restricted not to involve such lattice artifacts in order to keep proper connection with 
continuum theories. 

We have carried out a similar analysis in four dimensions by taking the link variables 
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(17) 



where s is an integer (|s| < jl^ corresponding to 77134 = s. Only the hnk variables along 
the first two lattice coordinates are changed. Since the spectrum of H is symmetric at t = 0, 
we can find indexD as in two dimensions. The characteristic features of the spectral flows 
observed in two dimensions can be seen also in four dimensions. In Figure |^ we indicate the 
spectrum of if for L = 4, |A| < 1, s = 1 and \t\ < 8. Most of the eigenvalues are distributed 
outside the parallelogram region and an eigenvalue crosses the parallelogram downward as 
t increases by unity from an integer. When it crosses the horizontal axis, indexD increases 
by one unit. Note that the direction of the eigenvalue flows on the parallelogram region is 
consistent with the factor (—1)^ on the rhs of (piSf) . The characteristic feature of the spectrum 
is not changed for s > 1. The parallelogram region becomes smaller as s increases. In general 
s adjacent eigenvalues flow downward and indexD increases by s when they cross the t-axis. 
This is consistent with (0). We found that the index and the topological charge coincide 
with each other for the values of (s, t) {0 < s < t) given in Table 0. Though our analysis is 
restricted to rather small lattice sizes 2 < L < 4, we anticipate that the parallelogram region 
expands rapidly enough as L increase and the condition (|^) with < 77 < (2 — -\/2)/12 is 
sufficient to ensure the lattice index theorem (|13]) in four dimensions.]^ 

^For gauge fields satisfying (||) and ( p^ nonvanishing topological charges can be realized only for L > 9. 
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none (1,1) (1,1), (1,2), (1,3), (1,4), (2,2), (2,3) 



Table 2: 



Coming back to two dimensions, we give some analytic results on the spectrum of H 
and indexZ) at t = r, L^/r, where r is an arbitrary integer factor of L. They are helpful to 
understand our numerical results. 

The eigenvalues A of if at t = L^/r with L a multiple of r are parametrized by p and q 



27ik 



2rnl 



P 



L 
2 



<k< 



L 
2 



< / < 



and satisfy the secular equation ^ 

det 



/5(p,g)-A Cip,q) 
[ C{p,qy -B{p,q)-X^ 



where B{p, q) and C(p, q) are r x r matrices defined by 



{B{p,q))ki 
{C{p,q))ki 



(a) , 

k,l+l "T 



2 . ^l-cos|^p + 
lsl^l, + sm(p+'^^ 



27Tk' 



(9) 

k+l,l 



(1^ 



(19) 



4^ + 2^tM- (0<fc,/<r-l) (20) 
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The 6^^] is the Kronecker's 5-symbol for < fc, / < r — 1 and satisfies the twisted boundary 
conditions S^^l = e~^'^6k,o and S^'^l = e'^'^o./t- It can be shown that (|T9|) takes the following 



form Il2i 



/.(A) = ^^sin^|sin^|, (21) 

where fr = A^*" + ■ ■ ■ is a polynomial of degree 2r and is independent of p and q. For 1 < r < 6 
it is given by 

/i = A^-l, 

/2 = X' -6X^ + 1, 

/4 = A^- 12A6 + 42A^-8A3-44A2 + 24A-3 , 



8 16 



(85 + 35V5)Vl0 + 2y5 109 + 50^5 

r A H — 



16 

h = A^2 - 18Ai° + ^A« - 6v^A^ - 360A6 + 54v^A5 

8145^, _ 297^3^3 _ 2007^, 255^3^ _ 719 _ 
16 2 8 2 16 ^ ' 

The case of r = 1 corresponds to the free spectrum since the link variables become trivial. 
The eigenvalues must satisfy the inequality 

< (-l)'-'2'^-Vr(A) < 1 . (23) 

This gives rise to 2r allowed intervals A^^-* < A < A^^-* (s = 1, ■ ■ ■ , 2r), where A^^-* and A^"^-* are 
the s-th largest roots of either /r(A) = or /r(A) = (— 1)''^^/2''~''. In general the allowed 
intervals become shorter and shorter as r increases. This is the reason why eigenvalues form 
clusters around the point where t = /r takes an integer. In particular we find (|TB]) from 
( p3|) for r = 2. Similar mechanism of clustering is also at work for other integer points t 
( 7^ 1) as is seen in Figure |l]. 

If fr is known, it is possible to find indexD at t = L^/r. As an example, we consider 
the case that L is a multiple of 4 and take r = 4. For any p, q the eight roots of ( pT]) are 
separately located in the following narrow intervals 

-2.4142 • ■ ■ < A < -2.4087 ■ ■ ■ , -1.9429 ■ ■ • < A < -1.9032 ■ ■ ■ , 

-1.7320 ■ ■ ■ < A < -1.6914 ■ ■ ■ , 0.1029 • ■ ■ < A < 0.1939 ■ ■ ■ , 

0.4142 ■ ■ ■ < A < 0.5609 ■ ■ ■ , 0.9260 ■ ■ ■ < A < 1 , 

1.7320 ■ ■ ■ < A < 1.7448 ■ ■ ■ , 2.7082 ■ ■ ■ < A < 2.7092 ■ • • . (24) 
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It has three negative roots and five positive roots for each p, g, giving the exact index 
— lxLxL/4 = — Note also that there is a large gap between the largest negative root 
and the smallest positive root. In general (pTD has r — 1 negative roots and r + 1 positive 
roots for r > 4 and, hence, indexD = — L^/r at t = L^/r. On the other hand it is vanishing 
for 1 < r < 3 since the number of positive roots and that of negative ones are always equal. 
These results can be extended to integers t = r = L/ s {< L) hy noting the relation 



det(/7-A)|i=, = (/,L(A)r , (25) 

where r and s are positive integers satisfying L = rs. This gives at t = r the index — r for 
sL > A and for sL < 3. These are completely consistent with the numerical results. 

We have confirmed that the equality (0) between the index and the topological charge 
holds true for a wide class of gauge fields on the finite periodic lattice in two and four 
dimensions. The condition (|^) with ([I2|) excludes uniformly the configurations for which 
the discrepancy between the index and the topological charge appears and ensures the in- 
dex theorem (^3|). We should give a remark on the possibility of extending our analysis 
to nonabelian theories with nontrivial topological charges. Link variables with nontrivial 
topological charges can be easily found within the commuting subgroup. The analysis can 
be reduced to that of a multi-U(l) cases and our results can be applied straightforwardly. 
The index theorem (|13]) also has an implication on chiral gauge theories on finite periodic 



lattices. We will argue them elsewhere. 

I would like to thank Peter Weisz for giving me the valuable suggestions to improve my 
primitive numerical results. I am also grateful to Dieter Maison and Theory Group for the 
kind hospitality during my stay at Max-Planck-Institute. 
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